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Abstract 

OO We consider a boundary-value problem for the second order elliptic differential opera- 

tor with rapidly oscillating coefficients in a domain Jig that is e— periodically perforated by 
psj small holes. The holes are divided into two e— periodical sets depending on the boundary 

^ interaction at their surfaces. Therefore, two different nonlinear Robin boundary condi- 

^ tions -|- eKm{us) = £gi"^\ m = 1,2, are given on the corresponding boundaries of 

the small holes. The asymptotic analysis of this problem is made as e ^ 0, namely the 
convergence theorem both for the solution and for the energy integral is proved without 
using extension operators, the asymptotic approximations both for the solution and for 
p ^ the energy integral are constructed and the corresponding error estimates are obtained. 

< 

1 Introduction and statement of the problem 

, In recent years, a rich collection of new results on asymptotic analysis of boundary- value 

problems in perforated domains is appeared (see for example [1]-[10]). The classical method 
^ proposed by E. Khruslov [H] and D. Cioranescu and J. Saint Jean Paulin [T2j is based on 

a special bounded extension of solutions in Sobolev spaces. It was established by V. Zhikov 
^ [S |9] that the homogenization results can be obtained without using the extension technique 

^ in Sobolev spaces in periodically perforated domains. It should be mentioned the paper |2], 

\Q where the homogenization results for an elliptic problem with a nonlinear boundary condition 

O in a perforated domain were obtained with the help of a new unfolding method that does not 

Q need any extension operators as well. 

In this paper we use this simple Zhikov's approach and the scheme of the paper [13], 
where the full asymptotic analysis (the convergence of the solution and the energy integral, 
the approximation for the solution and the corresponding asymptotic error estimate in the 
^ Sobolev space H^) was made for an elliptic problem with a nonlinear boundary condition in a 

thick junction. 

Let S be a finite union of smooth disjoint nontangent domains strictly lying in the unit 
square □ := G M" : < < 1, i = 1, n}. In an arbitrary way, we divide B into two sets, 

B^^'> = IJ B^^^ and B^'^^ = |J bI^\ Let us introduce the following notations: 

k=l k=l 

Qo:=D\ B, B^""^ := [j^^^^ {z + E^™)) , iS^™) := eB^"^^ = {xeW: e^'x G B^™)}, 

m = 1,2, 



X 
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where e is a small parameter. Let f2 be a smooth bounded domain in M". Define the following 



perforated domain fig = \ (^Be U Be ) and require the domain fi^ to be a domain with the 

Lipschitz boundary. Denote = dQ fl fie and E, 
(see Fig. 1). 



(m) 

■£ 



QndB. 



(m) 
e 1 



m = l,2, H, = H(^)uHf. 




Figure 1: 



Let aij{^), ^ G M", i,j = l,n, be smooth 1— periodic functions such that 

1) V2, J = 1, . . . , n, G : a,,(0 = ^,.(0, 
2) 3xi>0 3x2 >0 V^gM" Vt^gM": xilr^]^ < aij(0^i^j < ^2!^/^ (1) 

Remark 1. Here and in the sequel we adopt the Einstein convention of summation over 
repeated indexes. 



Let fe, fo,g£^\glr^ be given functions such that fe, fo G L'^{fl), ge"^' ,go^' G H^{Q) and 

fe^fo in L\n), gij-)^g{n.) ^^^^^^ ^i(^)^ m = 1, 2. (2) 



(m) (m) 



The given functions Hj. 



m 



Km G Wl^^iM.)) and such that 

3 ci > 3 C2 > : ci < k'^< C2 a.e. in M (m = 1, 2) 



1, 2, are Lipschitz continuous (it is equivalent that 

(3) 



In the perforated domain Vte we consider the following nonlinear problem 





= fe 


in 


a, 


ae{Ue) + eKi{Ue) 


= egi'' 


on 




aeiUe) +eK2iUe) 


= ec/P 


on 


■^(2) 


Us 


= 


on 


Te, 



(4) 



du 



where Ce{ue) = d^^ {alj{x)d^Me{x)) , (Te{ue) = a'lj{x)d^Me{x) Ui, afj{x) = aij (|) , d^M - q^^, 
(z/i(|), . . . '^n(f)) is the outward normal. 
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Recall that a function Us from the Sobolev space H^{Qe, ^e) = {u E H^{Qs) '■ u\r^ = 0} is 
a weak solution to problem Q if the following integral identity 



la — /^C"") /o ^ — ^ /^C"") 



(5) 



holds for any function ip G H^{Q£, 



Our goal is to study the asymptotic behavior of -u^ as e — 0. Also it will be understandable 
further how conduct research in the case of p-multiphase interactions in perforated domains. 



2 Auxiliary uniform estimates 



Let HUQo) = {ve H\Qo) 



V — 1-periodic in ^i, . . . ,^n}- Obviously, we can periodically 
extend every function v from H^^^iQa) into Hl^[R'' \ {BW Ui3(2))); this extension will be 

denoted again by v. Let ip^^ G H^^j.{Qq), m = 1,2, be weak solutions to the corresponding 
problems 



per 

= qi in Qo, 
■1 on 
:0 on 
0, 



' ^id'^o^) = 92 in Qo, 
(T^(42)) = on ^W, 

a,(^(^)) = l on ^(2), 



(6) 



where = % (0^^(055^^ (0) , o-e(V') = aijd^^ip{^)iyi{^), (z/i, . . . , i/„) is the outward 

normal to S = S^^^ U S^'^\ S^™) = = l^^"^)] = meass^^™) (m = 1,2), 

\Qo\ = meassQo, ('^)qo — iQo'^iO The existence and uniqueness of the solutions to 
problems (|6| follows from the lemma. 



Lemma 1. Let Fi G L\Qo),i = 0,n, F^^j G L2(^( 
G ifpgj,((5o) ^0 following problem 



m 



1, 2. There exists a solution 



a^(N) = -F,u, + Fi% 



on 
on 



^(2) 



(7) 



if and only if 



In addition this solution is defined up to an additive constant. 



The proof is standard (see for instance [5J). Then the e-periodic functions tp^^ (-) 
Qe, m = 1,2, satisfy the following relations 

(4^^ (f)) 



X G 



dx 



' (4 



u v^) axj 



(4'Hf)))=^-'?2 in a 



-1 "(1) 
e ^ on ', 



on si'), 
on ^£ ' . 



e 



'1 



Multiplying with arbitrary function Lp G if^(fie,re) the corresponding differential equation, 
integrating over VL^ and taking into account the boundary conditions, we get the following 
integral identities 

el ifdsx = e alj{x) di:.^^^\C)\^=?i d^^ip dx + qm i fdx m = 1,2. (9) 

Due to the regularity properties of solutions to elliptic problems we have 

sup|V54"^(Ok=f I = sup I Vg^"^ (01 < Co (m = 1,2). (10) 



Using Cauchy's inequality with S {ah < 5a^ + i^; a,b,5 > 0) and (10), we deduce from (9) the 
following estimates (m=l, 2) 



'pds^<Ci[£'' \Vcc'p\dx+ Lf'^dx], (11) 



j ^^dx<C2(^e^j dx + e j^^^ ds}j W^eH\Qe,'re 



(12) 



where the constant Ci and C2 are independent of e. 

Remark 2. In what follows all constants {Ci} and {q} in inequalities are independent of the 
parameter e. 

It follows from ^ and (g that 

tLii^^^il^(.<-))<^3. (13) 



Also with the help of (11 ) and (|12j) it is easy to prove that the usual norm || ■ is uniformly 

equivalent with respect to e to a new norm 

1/2 



■wile := ( / \Vu\^ dx + e I ds. 



He 



in the space H^{Qs, T^), i.e., there exist constants C3 > 0, C4 > and > such that for any 
e e (0,eo) and u G H^^Qsj^e) the following relations hold 

C3||w||Hi(f7,) < \\u\\s < C^WuWmin,)- (14) 



2.1 Existence and uniqueness of the solution to problem ([4|) 

Associated with (|4]), we consider the energy functional 

h[u] := \ j al^{x)d,^ud,^udx + eY,{ j K^'^\u)ds,- j g^^^^uds,) - j f^udx (15) 

He m—1 ^^^j --(m) He 

on H^{VL^,T^), where 

fs:(™)(^) = ! K^{t)dt V^gM, m = l,2. (16) 
Jo 

It is easy to prove that if is a minimizer of at a fixed value of e, then is a weak 
solution to problem 



Theorem 1. At each fixed value of e problem has exactly one solution Ug G H^{^li;,Ti;) fi 
which the following estimate 



or 



(17) 



holds, where the constants Ci and C2 are independent of e, f£,gi"^^ and Ug. 
Proof. Integrating inequalities in (|3|, we obtain 

Clt^ + f^miO) t < Kmit) t < C2t^ + Km(0) t Vt G M, 

whence it fohows that 

^z'^ + KraiO)z<K^"'\z)<^Z^ + KmiO)z \/ Z E R m=l,2. 



(19) 



Using (14), (18), (19), ^ and the same arguments as in Theorem 1 ([13]), we can prove 
the coercitivity condition on /, i.e., the following inequality 



(20) 



holds for any function u G H^{Qg, F^). 

With the help of (g we can re-write the energy functional as 



+ q^j K^"'\u)dx-e J <A(^o(0)k=fHx.^?f)+^?f)M dx- 



He 



^r'^udx 



firU dx. 



Consider the function 



Hp, = \(AjPiP3 + (^4- (^0(6) k=f i^m{t) Pi + qmK^'^\t)- 

- ea%di, (^o(O) k=f [td^S^r^ + - ^..tgf^^) - fet. 



m=l 



Since 



9p,p^L{P, ^) ViVj = 2 alj{x)r]ir]j > xi Vp, r] G M", x G n^, 

the function L is uniformly convex in p for each x G fi^. This means that /[■] is weakly lower 
semicontinuous on H^{Qg,rs) and there exists at least one minimizer (see [T^ Chapter 8.2]). 

Thanks to ^ it is easy to prove the uniqueness of this minimizer (see Theorem 1 (jl3j)). 

Finally, let us deduce the uniform estimate (17). Denote by Ug the solution to problem 
Setting (y9 = in (|5]) and taking into account ([ij) and the left inequality in (18), we get 

^1 / iVusl"^ dx + e ci / ulds^ + eKm{0)y^ / Usds^< / fsU^dx + ey^ / g'^^^u^ds^ 

In I" /=(»") In 



from which 



{Q.,)\We\\L'^{^^) + 



m=l 



Using (|14|) and (|ll|), we derive the first part of the estimate (17) from the last inequahty, and 

□ 



then the second one on the basis of ^ and (13). 



3 Convergence theorem 

In the sequel, y denotes the zero-extension of a function y defined on fi^ into the domain Vl. 
Also we introduce the following characteristic function 

[0, xeU\QQ. 

It is known that Xqo(^) XQo(f) IQo| weakly in LF'iVL) as e — > 0. 

Lemma 2. Let {t'e}e>o he a sequence in if^(fie,re) uniformly hounded in e in H^{VLi.,Vs) and 
such that 

i^m{v^) ^ C, weakly in L^{Vl) as e^O (m = l,2). 
Then for any function ip G H^{VLi;, F^) 

e/ K,m{ve)^dsr, ^ Qm ({x) ip{x) dx as E^O (m = l,2). (22) 

Proof. By virtue of ([9]) we have 



(m) 



£ / k{vs) ipds^=e / a-j-(x) 95^Y'o™^(Ok=f {'^'iPe) dx.Ve ip + k{vs) d^^p) dx+ 



+ qm nive) p dx, m = 1, 2. 



Thanks to the Lemma's condition, ^ and (10), the first summand vanishes and the second 
one tends to qm ({x) p dx as e m = 1,2. □ 



Remark 3. From Lemma^ it follows that for any sequence {f£}e>o G H^i^ei'^s)-, which is 
uniformly hounded with respect to e, there exists a suhsequence {e'} C {e} (again denoted by 



{e}) and a function ( G L'^{Q) such that the convergences (22) hold. 



Using (K, we can prove similarly as in Lemma [2] that for any function p G H^{^1^, F^) 

(""^pds, ^ \S^"''^\ [ g^o^\x)pix)dx as e^O (m = l,2). (23) 



-,(m) 
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Consider 1— periodic solutions T;, / = 1, . . . ,n, to the following problems 

cr^{Ti) = -auUi on 5", (T;)qo = 0. 

From Lemma [T] it follows the existence and uniqueness of the solutions to these problems. 

With the help of T/, / = 1, . . . , n, we define the coefficients of the homogenized matrix {aij} 
by the formula 

aij = {aij + aikd^J'j)Q^, i,j G {1,2, . . . ,n}. (25) 

It is easy to see that 

a^J = {an % (6 + T.) d^, (e, + Tj))Q, (26) 
i.e., the matrix {aij} is symmetric and it is well known that it is elliptic (see for instance ^). 

Theorem 2. For the solution Ue to problem there exists the following convergences 

Ue ——^ IQol'^o weakly in L'^{^), 
afjdxjUs ttijdx^Vo weakly in L'^{^), i = l,...,n, 

where vq is a unique weak solution to the following problem 

2 2 



as e ^0, (27) 



-a^J^l,^Vo{x)+ j:\S^"^^\Km{Vo{x)) = E IS^""^ \ 9t\^) + \Qo\ Mx) , X G 
m=l m=l 

Vo{x) = 0, X E dQ, 



which is called homogenized problem for Q). 

Furthermore, the following energy convergence holds as e ^ : 



(2^ 



/ aij d^^vo da:,vo dx + V'lS'^™^! / Km{vo) vq dx =: Eo{vo) . (29) 
Jn Jn 



Proof. 1. It follows from (17) and (iSh that the values 



ll«e||L2(n), \\aljd^^u^\\L2(n), i = l,...,n, \\Km{ue)\\L^n), m = l,2, 

are uniformly bounded with respect to e. Hence there exists a subsequence {e'} C {e}, again 
denoted by {e}, such that 



Me IQol'^o weakly in L^(fi), 

aij d^Me 7i weakly in L'^{Q), i = l,...,n 

f^miue) — Cm Weakly in L'^{^), m = l,2. 



> as £ ^ 0, (30) 



where Vq, ji, i = 1, . . . ,n, Cm-, rn = 1,2., are some functions which will be determined in what 
follows. 
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2. Obviously the e-periodic functions T; (-) , / = l,...,n, defined in (24) satisfy the 
following relations 

d^. K(e)55,TKOk=f ) + d^M^) = Vx G a, 
Multiplying the first relation by (f), where (f) is arbitrary function from C^{Q), and integrating 



over fig, we obtain 



/ {ciijiO di^jTiiO + (^iiiO)k=f {ued:,,^(f) + (pdx.Ui;) dx = 0, 1 = 1, n. 



(31) 



Put the following test-function = eT;(|)0(x), x G fi^, into the integral identity ([5j). 
The result is as follows 



(^lji.x)dx,Ued£^Ti{C}k='^ <p{.x) dx + e alj{x)dx^Ue Ti{-) d^^(f){x) dx+ 

+ V / K^{ue)Ti(t>ds^ = e feTi<f)dx + e^S2 g^r'^Ti<pds,. (32) 
Using (|2|), (g and the identities (|9|, it follows from (g that 



/ afj(x)9:i..M£%Ti(^)|g=| 0(x) dx = 0{e) as e ^ 0, I = l,n. 



(33) 



Subtracting (32) from (31), we get 



{aij{^)d^Ti{^) + aii{^))\^=!LUed^^(j)dx + / alid^^.Ue (p dx = 0(e), l = l,n. (34) 



In (34) we regard that the functions aij d^.Ti + an, I = 1, . . . ,n, are equal to zero on B. 



Let us find the limit of the first summand in the left-hand side of (34). At first we note 



that the limit function vq in (30) belongs to HliVt) because of the conectedness of the domain 



W \ (i3(i) US(2)) (see [8]-[l0]). Since (aij(0%Tz(0 + aii{^))vi{^) = at ^ G 5 and the 
vector-functions 



(35) 



are solenoidal in Qq (see (24)), their zero-extensions into □ \ Qo are also solenoidal in weak 
sense, i.e., 

/ Fi{0-vmd^= [ Fi{0-Vij{Od^ = o ^^ec^jn), l = l,...,n. 

Then using results by V.V. Zhikov (see [H', Th. 2.1]), we get that 

lira / {aij{^)di:Ti{^) + aii{^))\i:=!iUsd:^^(f)dx = / au Vo dx^(f) dx . 



As a results, it follows from (34) in the limit passage as e ^ that 



aiiVod^^(f)dx + / -fi(f)dx = V0 G C^(fi), (/ = !,. ..,n), 



I.e., 



'yi{x) = an dx^vo{x) for a.e. x e (/ = l,...,n). (36) 

4. Using the extension by zero and the identities (|9]), we rewrite the integral identity ^ 
in the following way 



alj dxjUs d:r,f dx+ 



m=l 



+ qmjj_ho9^:^^^d^ V^gCo-(1]). (37) 

It is easy to see that the pointed summands in (37) vanish as e — > 0; the first one due to 
(|3|, (|T0|) and (|l7l), the second one due to (|To|) and ^ 

Taking into account (30), (36) and (|2]), we pass to the limit in (37) as e ^ 0. As a result 
we get the identity 

/ aijdx^vod^^Lpdx +^^qm I Cm^dx =\Qo\ / fo^p dx + ^^\S^'^'^\ / gl^^^ipdx (38) 
Jn Jn Jn Jn 



for any function ip G C^{Q). Since the space C^{Q) is dense in Hq^Q), identity (38) is valid 
for any function (f G Hq{Q). 

5. With the help of (|2|), (|5|) and (|38|) we can find that 



lim £*£(«£) = lim / feUedx + S^ j ^f^™^ rfs^^. = lim / feUedx+ 

+ Ef^ / (^U'')diA"'^k=^- {d^dr~^u, + gf^d,,M,)dx + q^ [ gt^^u^dx 



\Qo\ / fovodx + y^\S^'"^''\ I g^"'\odx= / a^dx^vodx^vodx + y^qm / 



(x)vq dx. 



(39) 



6. Now it remains to determine the last summand in (39). For this we will use the method 



of Browder and Minty, a remarkable technique which somehow applies to the corresponding 
inequality of monotonicity to justify passing to a weak limit within a nonlinearity. 
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Thanks to ([T]) and ([s]), the inequahty of monotonicity in our case reads as follows 
/ a^j d^^ {us-(p~ eTp d^^ip) d^^ {u^ - (p - eTg d^^ip) dx+ 

+ sJ2 [ ^^^{'^rniUe)-f^miv)){Ue-ip)ds,>0 \/ipeC^in), (40) 

m=l "^^^e"" 

which is equivalent to 



+ / atj {dxj<^ + d^Tp d^^ip) {d^^p) + d^Tq d^^ip) dx - 2 a^- d^M^ d^^p^dx- 
Jfie Jn 

-2 a\j d^^Ue d(.Tg d^^p^ dx~2e I a^j {d^^u^ - d^^p) - %Tp d^^cp) dl^^^p:>dx+ 

+ alj Tp Tg d^^^^ip d^^^p^dx- 

2 r 

(fim{p>) Ue + KmiUe) ip - Km{p>) <P>) dSx > Wp>eC^{^l). (41) 

m=l "'"E 



The limit of the first line in (41) is equal to the right-hand side in (39). The first integral in 



the second line can be re-written in the form 



/ (^P + + ^g)) k=f dx^^ dx- 



(42) 



It follows from [S] that its limit equals J^apg d^pP^ dx^p^ dx. Due to (33) the integral in third line 
vanishes. Obviously, the limits of summands in the fourth line are equal to zero. The limits of 
the integrals in the last line can be found with the help of Lemma [2] As a results we have 

/ aij dx^{vo - p>) dx^vo - </?) / (Cm - \Qo\K.7n{p^)) {vo - p^) dx > 0. (43) 

Jn Jn 

Evidently, this inequality holds for any function p) G Hq{Q). 



Fix any ijj G C^{^) and set p? := Vo - Xip (A > 0) in (43). We get then 



A / aij dxjip dx^t/j dx + y^^lm / (Cm - \Qo\K.m{vo - Xip)) ijjdx>0 yipeC^iO.). 
Jn Jn 



In the limit (as A — > 0) we obtain 



/ 9m (Cm - \Qo\KmiVo)) dx > 0. 
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Replacing ip by —ip, we deduce that in fact quality holds above. Thus 



QmCmix) = ^ l^^"^! K„Xvoix)) for a.c. X en. 



(44) 



m=l 



m=l 



7. Returning to (38), we see that the function vq satisfies the following integral identity 

2 



Jn Jn Jn Jn 



(45) 



for any function (p G Hq{Q). Hence Vq is a weak solution to the limit problem (28). Thanks to 
(|3| this solution is unique. 



Due to the uniqueness of the solution to problem (28), the above argumentations hold for 



any subsequence of {e} chosen at the beginning of the proof. By replacing (44) in (39), one 
obtains the convergence of energies (29). □ 



4 Asymptotic approximation to the solution 
and the energy integral 

We take the following approximation 

Ue := Vo{x) + eTk (f ) ^^^^^0(3;) 



(46) 



to the solution u^- Substituting the difference — u^, we find the residuals both in the differ- 
ential equation and boundary conditions. Straightforward calculation show that 



-Ce {Ue - Ue) = fe{x) - fo{x) - ^ Qm (5'o"'^(a;) - Kni{vo{x))) + 

m=l 

+ (aij(0 + aifc(0%^i(0 - m^ij) k=f ^Ix.vo + ed^^ {F[{x)) , x e Qe 

ae{ue-Us) = -eKm{ue)+egi"'\x)-F^'{x)iyi, xeS(™) (m = l,2), 
where Ff{x) = aij(f)Tfc(f ) d^^^^Voix), i = 1, . . . , n, and 



\Ur_ - U 



(47) 
(48) 

(49) 



Let be a smooth function in VL such that 1, v^e(x) = 1 if dist{x,dQ) < e, 

and (pe{x) = if dist{x,dfl) > 2e. Obviously, 



|Vx</2e| < ce ^ in Q. 
With the help of ip^ we define the following functions 

ipeix) = -eips{x) Tfc (1) d^^Vo{x) and We{x) = u^{x) - Ue{x) - tpeix), X eQe 



(50) 
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It is easy to verify that supp {ipe) C = {x G fi^ : dist(x, dQ) < 2e} and We is a solution to 
the following problem 

+ (ay(0 + a,fc(^)%Tj(0 - \Qo\~^aij) d^^^Vo in Q,; 

a,{ws) = -eK^{ue) + egi"^\x) - F[{x)ui- asitpe) on H^"^ (m = l,2); 
We = on r^. 

Multiplying the equation of this problem by w^, then integrating by parts and subtracting 
identities ^ for Lfm = Umivo) Ws, m = 1,2, we get 

m=l "^--^^ -^^^ 

= {fs - fo)<^dx + y^(e g^""^ Weds^ - qm g^^^Wedx)- 

2 r 

- £ ^ / a\^{x) 95^.^^015=1 ^x, {l^m{v^) We) dx+ 
m=l 

+ / i^ijiO + «ifc(0%^i(0 - \Qo\~%j) |e=f dlx,voWedx+ 

+ 6 F-d:r,Wedx~ / a^, d^^^pe dx.We dx . (51) 



Due to Q, ([3]) and (14) the left-hand side of (51) is estimated by the following way 

/ d^^We d^^We dx + e^ / {UmiUs) - 'tm(fo)) We ds^ > 

\Vwe\^dx + e w'^ds^^-C2e \{eTkdx^vo + iJe)we\dsx > 

> csllweWm^n^) - / \{eTkd^^^vo + iJe)we\ ds^. (52) 

Now estimate the summands in the right-hand side of (51 ). Evidently, | {fe ~ fo) Wedx\ < 
life — fo\\L'2{ne)\\we\\H^ne)- With the help of ([9|, ^ and (10) we bound the second and third 
terms: 



m=l 



^ / g^""^ We ds^ -qm fl-o™^ We dx^ 

/ gi^'^We dx- g^^^We dx) = Cie\\we\\min,) + C2 llfi-l 



£1 / 4- 9ijM0\.=^ d^X9i""^We) dx 



+ 



E 



m=l 



al-{x) di.^ipo\^^2: dx, inm{vo)we) dx 



< 
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< ecs / \Vvo\\we\dx + ec^'^ / \Kmivo)\\\/we\dx < ec^WweWm 



Thank to (25) and the fact that the vector-functions (35) are weak solenoidal in □, it follows 



from Lemma 16.4 ([3j) that 



in,)- 



It is easy to see that £\ F^dx^Wedxl < £CQ\\we\\m(n^). The last summand in (51) is 
estimated with the help of Lemma 1.5 ([S]) and (50): 



aljdxj^edx.We dx 



/ aljd^^tpsdx.Wedx 

JU2e 



< 



< Cy |Vfo| iVifel dx + e^cs / l-D^fol iVtWel dx < 
< C7\\vQ\\m(pi^^)\\we\\m(}A2e) + ^'^cs\\vQ\\H'i{n)\\we\\m{n^) < cge'^ \\vo\\H^(^n)\\w£\\Hi(^Q^). 



It is remain to bound the last term in (52). Thanks to ^ and (10) we have 



dx+ 



(53) 



+ cioe 



Tkd^.VQ We 



dx < ciie\\vQ\\H2in)\\we\\m(n,)- 



With the same arguments as in ( 53 ) and ( 54 ) we have 



e I \i'eWe\dS:r <Ci2e-^\\vo\\H2{n)\\We\\m{U,)- 



Finally, we conclude from (51), (52) and estimates obtained above that 

2 



\W^ 



m=l 



Since ||V'e||//i(f^£) bounded above by €26^, we have from (55) that 



(54) 



(55) 



\Ue — Up 



\min.) < C,{e'^ + \\fe- foh^in.) + H^^^ " dt^Wi^.)]. 



m=l 



where the constant C3 is independent of e. 
Thus, we have proved the following result. 



(56) 



Theorem 3. Between the solution Ue to problem ^ and the approximation function (4-6) the 
estimate (56) holds. 



With the help of the approximation function (46) and estimate (56) we can obtain an 
estimate for the energy integrals. 
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Corollary 1. The following estimate 



(57) 



m=l 



is satisfied, where the energy integrals Es{ue) and Eq{vq) are defined in (29). 



Proof. By virtue of ( 56 ) we have 
where 

2 

||rniL^(n.) < C^{e'^ + II/. - /o||l^(o.) + J2 H^^^ " dlT^hnn.) 



m=l 



Then 



(5^ 



where 



Pe = 2 a?, dx^vo + d^TkiO L=i d^^Vo )rldx+ / af,T*r^ dx. 



Taking into account the boundedness of afj and d^.Ti{^) and estimate (56), we get 

|Pe| < Ci I ||fo||Hi(n) 



r*r* dx \ + rlrl dx \ < 



< C2 + ll/e - foWmn^) + U""^ - air^hnn^) 



(59) 



Due to (2) we can regard here that ||/e"~/o||^2(Q^-) ^ ||/e — /o||L2(f2£)) similar for other summands. 
Let us introduce the following functions 

^fcKO = «..(e)%(Tfc(0+6)%(TKO+eO - ^flfc^, k,l = l,...,n. 

After extending the functions aij, Tk, d^T^, k = l,...,n, by zero to nyQo, the functions 
Hki, k,l = 1, . . . ,n, will be 1-periodic with zero average over □. 



By the same way as we rewrote a summand in (41) (see (42)) and using the functions 



Hf^i, k,l = 1, . . . ,n, and (58), we obtain 



/ aljdx^e dxMe dx - aijdx^vo d^^Vo dx = HkiiO l^=£ d^j^Vo dx^Vo dx+ 
JQ.^ Jn Jn " 

+ / ( i^XQod ) - l)aijdx^vo dx.vo dx + Pe =■ h + h + Pe, (60) 



where XQo the characteristic function defined in (|21j). The summand Ji can be estimated 
by the same way as in the proof of Theorem 1.3 ([51 Ch. 2]). As a result, we have |/i| < 

ci£|ko||^2(n)- 
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To estimate I2 we note that /□(|^XQo(0 " l)*^^ = 0- Therefore, with the help of Lem- 



ma 1.1 ([3j) we get I/2I < C2£:||t;o||j^-2 



Summarizing (59) and estimates for Ii and I2, from (60) we deduce as follows 



Jn^ Jn 



< 



(61) 



m=l 



Now consider the difference 

2 



With the help of integral identities (|9| we re-write it in the form 

h = ey^j alj d^^^''^\C) d^^{km{ue)ue) dx+ 

+ ^qm i km{Us)Ue - ^ |5""| / k^Mvo dx. 
m=l -^^^ m=l 

Due ([3]), (10) and (17) the first term is not grater then eci. Since 



krr,{u^)Ue dx - / k^{vo)Vo dx 



it remains to estimate the following difference 
2 



< C2\\u^ - ^^o||l2(0,), 



km / k„i{vo)Vodx - \S"^\ / km{vo)vo dx\ = {(ImXQoi^) - \S"^\)^rn{Vo)Vodx 

Thanks to the equality /□(o'mXQolO ~ \S'^'\)di = {qm = \S'^\/\Qq\) and Lemma 1.1 ([3j), 
this difference is bounded by C3e||wo||//i(c)- Thus, I/3I < c^^e + C2\\ue — t'o||L2(Qs)- 



Finally, taking into account the previous estimate, (61), (56) and noting that Ee{u, 



Eo{vo) = h + I2 + h + Pe, we arrive to (57). 



□ 
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